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Abstract 

A well-known method of transferring the population of a quantum system from an eigenspace of the free Hamiltonian to another 
is to use a periodic control law with an angular frequency equal to the difference of the eigenvalues. For finite dimensional 
quantum systems, the classical theory of averaging provides a rigorous explanation of this experimentally validated result. 
This paper extends this finite dimensional result, known as the Rotating Wave Approximation, to infinite dimensional systems 
and provides explicit convergence estimates. 
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1 Introduction 

1.1 Effective control of quantum systems 

The state of a quantum system evolving on a finite di- 
mensional Riemannian manifold f2 is described by its 
wave function, that is, a point in the unit sphere of 
L 2 (£l, C). In the absence of interaction with the environ- 
ment and with a suitable choice of units, the time evo- 
lution of the wave function is given by the Schrodinger 
equation 



.dip 

l at 



i 



Aip + V(x)ip(x,t), 



where A is the Laplace-Beltrami operator on ft (with 
suitable boundary conditions) and V : f2 — > R is a real 
function (usually called potential) accounting for the 
physical properties of the system. When subjected to an 
excitation by an external electric field (e.g. a laser), the 
Schrodinger equation reads 



M i A/ 

l— — = — Aw 

dt 2 v 



V(x)ij(x, t) + u(t)W{x)i>{x, t), (1) 



where W : fi —> R is a real function accounting for the 
physical properties of the laser and u is a real function 
of the time accounting for the intensity of the laser. 
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A natural question, with many practical implications, is 
whether there exists a control u that steers the quantum 
system from a given initial wave function to a given tar- 
get wave function (controllability issue) and, more im- 
portant, how to build this control law (effective design 
of controls) . 

Considerable effort has been expended by different com- 
munities on studying the controllability of (1). We refer 
to Nersessyan [NerlO], Beauchard & Mirrahimi [BM09], 
Mirrahimi [Mir09], Boscain & Laurent [BL10] and 
Boscain, Caponigro, Chambrion & Sigalotti [BCCS12] 
for a description of the known theoretical results con- 
cerning the existence of controls steering a given source 
to a given target. As proved by Nersessyan [NerlO] and 
Mason & Sigalotti [MS10], approximate controllability 
is a generic property for systems of the type (1). 

A number of effective control algorithms have been ob- 
tained by various authors, see among many others War- 
ren, Rabitz & Dahleh [WRD93], Bliimel, Fishman k 
Smilansky [BFS86] , Ohtsuki, Kono & Fujiyama [OKF98] 
or Belhadj, Salomon & Turinici [BST08]. Most of the 
controls used in practice exhibit a remarkable pattern of 
periodic shape, with a frequency corresponding to the 
transitions of the quantum system (see for instance Sa- 
lomon, Dion & Turinici [SDT05]). 

1.2 Averaging techniques in quantum mechanics 

The fact that a small amplitude periodic excitation with 
suitable frequency is sufficient, in general, to induce a 
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transfer of quantum population from one energy level to 
another has prompted much attention. 

The situation is now well understood for quantum sys- 
tems with finitely many energy levels. These systems ap- 
pear, for instance, as truncations of infinite dimensional 
systems. In this case, system (1) reads x = (A + uB)x 
where A and B are N x N skew-hermitian matrices. Let 
us briefly recall the method of the proof. 

The mathematical concept of averaging of dynami- 
cal systems was introduced more than a century ago 
and has now developed into a well-established the- 
ory, see for instance the books of Guckenheimer & 
Holmes [GH83], Bullo & Lewis [BL05] or Sanders, Ver- 
hulst & Murdock [SVM07] . It was observed that, for 
regular F and small e, the trajectories of the system 
x = sF(x,t,e) remain e close, for time of order 1/e, to 
the trajectories of the average system x = F(x) where 
F{x) = lim^ 1/t J*F{x,t,0). 

In quantum physics, this procedure has been known as 
the rotating wave approximation for several decades, see 
Fox & Eidson [FE87] or Vandersypen & Chuang [VC05]. 
Assuming without loss of generality that A is diagonal 
with eigenvalues iAi, . . . , iXn, define, for every n in N, 
y n : t i— > e~ tA x n where x n is the solution of x = (A + 
u*(t)B/n)x. The mapping y n is absolutely continuous 
and satisfies y n = (u*(t)/n)e~ At Be At y n . The conclusion 
follows from standard averaging theory by computing 
the average matrix lim^oo l/(nt) f Q u*(T)e~ Al ' Be Al "dr 
whose all entries are zero but maybe entry (j, k) if u* is 
2n/\\j — Afe|-pcriodic. 

By contrast with this finite dimensional result, the situ- 
ation is much more intricate when the ambiant space has 
infinite dimension and the problem involves unbounded 
operators, which is precisely the case for the Schrodinger 
equation (1) when dim 51 > 1 and A = iA/2. The above 
averaging method fails because of serious regularity is- 
sues. For instance, the mapping t H» e tA is no longer 
Lipschitz continuous. If moreover B is unbounded, that 
is, if function W in (1) is not in L°°(fl, R), then, for 
any given t, the mapping x H> e~ tA Be tA x is not even 
continuous: this is what prevents the direct application 
of the averaging results in Banach spaces presented by 
Artstein [ArtlO] . For these reasons, most of the available 
averaging results for infinite dimensional quantum sys- 
tems deal with constant controls only, as in the papers 
of Kummcr [Kum71] or Scherer [Sch94]. 



1.3 Framework and notation 

We first reformulate the problem (1) in a more abstract 
framework. In a separable Hilbert space H endowed with 
norm || • || and Hilbert product (•,•}, we consider the 



following evolution problem: 

dtp 



dt 



= (A + u{t)B)ijj(t) u(t) e U 



(2) 



where (A, B, U) satisfies Assumption 1. 



Assumption 1 (A, B, U) is a triple where (A, B) is a 
pair of linear operators and U is a subset o/R such that 

(1) for every n in N, U C nil ; 

(2) A is skew-adjoint with domain D(A); 

(3) there exists an Hilbert basis (</>fc)fceN of H and a 
family (iXk)keN in iR such that A<fik — i\k4>k for 
every k inN; 

(4) B is skew-symmetric, possibly unbounded with do- 
main D(B); 

(5) for every k in N, 4>k belongs to D(B); 

(6) for every u in U, A+uB is essentially skew-adjoint. 

Assumption 1.6 ensures that, for every constant u in U, 
A+uB generates a group of unitary propagators. Hence, 
for every initial time to in R and initial condition ipo in 

H, for every piecewise constant control u taking value 
in U, we can define the solution t i-> T u (t, to)ipo of (2) 
taking value f/'o at time to- We simply note T u (t, £0)^0 = 
T^ when t = 0. 

Remark 1 To the best of the author's knowledge, in the 
general frame of Assumption 1, no definition of solu- 
tions of (2) is available for controls that are not piece- 
wise constant. With some extra regularity assumptions 
(for instance, when B is bounded, see Section 3.4), one 
can define the solution of (2) for more general controls. 

I. 4 Main result 

Definition 1 Let (A, B, U) satisfy Assumption 1. A 
point (j, k) of N 2 is a non-degenerate transition of 

(A, B) if (i) j ^ k, (ii) (cpj,B(j)k) ^ and (Hi) for every 
I, m in N, \Xj — Xk\ = \Xt — X m \ implies {j, k} — {I, m} 
or{j,k}n{l,m} = <D or{(j) U B4, m ) = 0. 

Theorem 1 Let (A, B, U) satisfy Assumption 1, (j, k) 
a non- degenerate transition of (A, B) and u* : R+ — > 
U be a piecewise constant function, periodic with period 



T = 



|A;i-A fc | 



. At 



that 



f (r)e i(A! - A ™ )r dr = 0/o 



every I, m such that |A; — X rn \ € (N \ {l})|Aj — A^l and 
{(t> U B(j) m ) ^ and {I, m} D {j, k} ^ 0. 

f T 

If / M*(r)e l(A ^ Afc)r dT ^ 0, then there exists T* > 
Jo 



such that the sequence 
to I as n tends to infinity. 



tends 



2 



1 . 5 Content of this paper 

This paper comprises three parts. The first one (Section 
2) is concerned with a finite dimensional version of The- 
orem 1. As already mentioned, the convergence result is 
classical. A time reparametrization (Section 2.2) is in- 
troduced, which leads to explicit estimates better than 
those currently available in the literature. 

The second part (Section 3) contains a general proof of 
Theorem 1, valid in the general framework of Assump- 
tion 1. It provides the first available time estimates for 
the approximate controllability of general bilinear quan- 
tum systems when the free Hamiltonian admits a dense 
family of eigenvectors. 

Finally, Section 4 presents an extension of Theorem 1 
to the case of the Morse quantum oscillator. The main 
feature of this system is that the spectrum of the free 
Hamiltonian has a continuous part. 



2 Finite dimensional estimates 

2.1 Finite dimensional framework 

Let N be an integer and A^ N \ B^ N ^ be two skew- 
adjoint matrices of order N . Without lost of generality, 
we may assume that A^ N ^ is diagonal with eigenval- 
ues iAi,iA 2 , . . . ,iAjv. We denote with (0j)i<j<Af the 
canonical basis of C N , with (-,-,) the canonical Hermi- 
tian product of C^, with (b jk = (4> 3 ■, 5 (JV) 0fc)) x < j k < N 



the entries of and with 



A*) 



iff 



C N the 



(AT) 



projection it 
components in C 



^ EnKll^. 1 )^ 011 the nrSt 1 



N 



For every piecewise constant function u : R — > R, we 
denote with X™ N ^ (t, s) the propagator between times s 
and t of the system 



x = {A {N) + u{t)B {N) )x{t). 



(3) 



Assume that u* is periodic with period T 



2tt 



|A 2 -Ai| 



and 



that J ii*(r)e i(Al A "* )T dr = for every {I, m) such that 

{J, m}°n {1,2)^0 and A; - A m e (Z \ {±l})(Ai - A 2 ) 
and k m ^0. 



If i u*(T)e^ x ' 2 Al - )T dr ^ 0, then, for every n mN, there 
Jo 

exists T* in (nT* - T, nT* + T) such that 

i - 1(02,^(^,0)01)1 (i + 011^^) || 



{l + 2K\\B( N )\\)I 



n 



with 



TFT 



2|6i )2 | J V(T)e i ( A i-^) r dT 



= f\u{ 
Jo 



r)|dr, 



IT* 
K=—,C 



sup 

Q\fc)eA 



J T M*(T)e i ( A J- A ^ r dr 



. |Aj— A fc | 
|A 2 -Ai| 



where A is the set of all pairs (j, k) in {1, ... , N} 2 such 
that bjk ^ and {j,k} n {1,2} and | Aj - A fe | <£ 
Z|A 2 -Ai|. 

Corollary 3 With the notations of Proposition 2, 
| (02, X^(nT*, O)0i ) | tends to 1 as n tends to infinity. 

The convergence result (Corollary 3) is classical. The 
novel element here is that the bilinear structure allows 
us to give explicit estimates (4) for the convergence rate. 

Most of the rest of this section is devoted to the proof 
of Proposition 2. A technical time reparametrization 
(Section 2.2) is introduced. The main ingredient of the 
proof is an averaging procedure (Section 2.3), that is 
performed first on piecewise constant controls and then 
extended to irregular controls. Finally, we introduce the 
notion of efficiency in Section 2.4. 



In other words, t n> X? N , (t, s)xq is the unique solution 

of (3) with initial condition x(s) = xq. It is classical that, 
for every t, s in R, for every x in C N , the mapping u i-> 
X^ N At, s)x admits a Lipschitz-continuous continuation 

to Lj oc (R, R). 

The aim of this section is to prove the following result. 

Proposition 2 Let u* : R + ->Riea locally integrable 
function. Assume that \\ ^ A 2 and that, for every l,m < 
N, | A; — A m | = |A 2 — Ai| implies {l 7 m} = {1,2} or 
b lm = or {l,m}n {1,2} = Q. 



2.2 Time reparametrization 

We note PC the set of the piecewise constant func- 
tions for which there exist two sequences {uj)\<j< p and 
(tj)i<j< p with value in (0, +oo) such that 

U = U jX[T„T,+t,), 
l<j<p+l 

where \ is the characteristic function and the sequence 
( T j)i<j<p+i i s defined by induction: t\ — 0, Tj + i = Tj + 
tj. An element u of PC will be denoted {uj,tj)i<j< p . 
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The involutary mapping V : (uj,tj)i<j< p G PC i-> 

(l/uj,Ujtj) 1< . <p e PC can be used as a time 

reparametrization to replace the system (3) by the 
control system in C N 



dx 
dt 



= («(f)#)+#))i(t), 



(5) 



whose propagator between time s and t will be denoted 
with X? N) {t,s). 



Finally, for every t, we define z n (t) — e~ Vn ^ A y n (t) 
(note that, for every t, for every I in N, \{<f>i, z n (t))\ = 
\(</>i,y n (t))\), and the time varying N x N matrix M n 



M n :t^ sgK o v n )e-^ A(N) B^e^ A(N \ 



From (7), we deduce the dynamics of z n , valid for almost 
every t in G+ U G~ : 



Proposition 4 For every xo in C and every u in PC, 

Proof: This follows from the equality e t (^ 4<JY) +« B(lv) ) — 

e v « + 7 , valid on every interval where u is con- 
stant. □ 



dZr, 
dt 



= M n {t)z n {t). 



(8) 



We note the propagator associated with (8). Note 
that, for every k, the mapping t i-> (<f>k,Zn) is Lipschitz 
continuous with Lipschitz constant || || . 



#.3 Averaging procedure 

Let u* be a non vanishing piecewise constant function, 
T-periodic, as in the hypotheses of Proposition 2. For 
every n in N, define the non- vanishing T-periodic func- 
tion u n = u* In. For every n in N, J Q ™ |u„(r)|dT = 
/ T |«*W|dr. 

We perform now the time reparametrization set out 
in Section 2.2. The function t n> J* \u n (r)\dT is non- 
decreasing. We denote with v n its reciprocal function 
and define also I = J^\u*(T)\dT. For every t in R + , 
v n (t + I In) = v n (t) + T and V\u n \ is the I/n periodic 
derivative (defined almost everywhere) of v n . 

Fix Xq in and note with t n> x n (t) the solution of 



x = (A^ + u n {t)B^)x(t) 
with initial condition x(0) = x . 



(6) 



The set [0, T\ can be written as a finite union of disjoint 
intervals 

[0, T] = U Jf U . . . U J+ U Jp, 

such that u* takes positive (respectively, negative) val- 
ues on J+ = Uf =1 J ; + (respectively, J~ = Uf =1 JT). 
Defining the sets G+ = v [ n 1] (J+) = {I e K+\3s e 
J + ,Jo\u n (T)\dT = 1} and G~ = v [ ~ 1] (J-) = {I e 
R + |3s e J t Jo \un(r)\dT = /}, we obtain the dynamics 
of y n =i„oD ni valid for almost every t: 



Let M* be the constant N x N matrix whose entries, for 
^ < 3,k < N, are defined by 



if T(Xj - A fe ) e 2ttZ and mt fe = if T(Xj - A fe ) £ 
27rZ. Define the N x N matrix H n (t), with entries 
(^ k (t))i<j,k<N, by ff n (t) = Af n (i) - Aft. 



Lemma 5 For every t m R+ 7 /or every n m N 7 for every 
j,k < N such that \Xj - X k \ Z|A 2 - Ai|, 



/' 

Jo 



h%{r)dr 



< 



n 



J T u*(r)e iT ^- x ^dT 



sin 



/ . 



(9) 



For every t in R + , /or every n in N, /or every j,k < N 
such that{j,k} ^ {1,2} and |Aj — Afc | / 1 A2 — Ai| € Z, 



^(r)dr 



< 



21 b jk \I 



(10) 



Proof: For every t, define the integer s = . For every 
j, k, n, t, we have 



dyn = { (V(\u n \)(t)A^ +B( N ))y n (t) if t G G+ 
d * I (P(|w„|)(t)AW -#>)!/„(() if t 6 G- 



(7) 



< 



4 



For every j, k such that T(Xj — A&) £ 27rZ, 

\ sgiunov^e^-^^dr 
Jo 

s r I/ n 

pi /n s 

/ sgKo^je^^)"^ 1 -) V e i(A ^ Afc)mT dr 

^ m=l 

/ V(T)e i (^- A '=) r dT , ' 



< 



< 



1 

< - 



n \ |l-exp(iT(Aj-A fc ))| 



for every n in N, for every t < Vn ^ (K), for every xq in 
C N , (I, K and C are denned as in Proposition 2) 



/(C+l)||Mt|| 
1 + A-(||Mt|| + su Ps<4 0)x \ 



< 



or 



(12) 



|Xfo(t,0) -e^ {N) e^ 1] ^\\ < 1 + 2 K\\B^ 



J(C+ 1)||BW|| 



'• (13) 



If T(Xj - A fe ) e 2?rZ, then 



V / sgKo Un )e i (^-^)^^)+ roT )dr 

fl/n 

= s sg(u n ov n )e i( - X *- X ^ v "^dT 
Jo 

Jo 



s 
n 



Hence 

si /n 



psl jn 

/ sgKo^y^-^-^dT 

Jo 

- - ( T u*{ T )e l ^- x ^ T &T 
1 Jo 



< -■ □ 

n 



Recall that, for every n in N and every t > 0, 

i„ = M f 2:„ + H n (t)z n , 
or, using the variation of the constant formula, 



z n (t)=e tM 'z n (0) + I 
Jo 



= (t-s)M t 



H n (s)z n (s)ds. 



Integrating by part, we get for every t > 0, 



z„(t) = e* Mt z„(0) + 



D (t-s)M f 



f H n (r)dT) z n (s) 

+ |Mt e < ! - S ) Mt ^ H n (T)dT^j z n ( S )ds 

- jf' e ( *- s)Mt (jf ff„(r)dr) z n (s)ds. (11) 

The equality ^(4,0) - e M<W) o Z^vl^it), 0) gives 
estimates for the convergence of the propagators X u ": 



Projecting (11) on the two first components of C N , we 
get 



4 N) z n (t)-4 N) e tM 'z n (0) 



< 



4 N) /V(T)dT 

Jo 

4 N) /*ff„(r)d 
Jo 



4 W) / ff„(r)dr 



(14) 



Note that \(<j> 2 , e KM Vi}| = 1- Equation (4) follows from 
(14) with t = K and 



IS =«„(#) = «* 



nirl 



2|6i, 2 | J ti*(r)e i ( A2_Al )TdT 



Knowing that w* is non-decreasing and v*(ll) = IT for 
every Z in N, we deduce that, for every n in N, nT* — T< 
T* < nT* + T, with T* defined as in the statement of 
Proposition 2. 

Note finally that, for every s, t in R such that s < t, 

\\B<h\ 



\\{<h,X n (t))\-\(<h,X n (8))\\ < 



\u(r)dT. 



This completes the proofs of Proposition 2 and Corol- 
lary 3 in the case where u* is a non vanishing piecewise 
constant function. 

If u* is a locally integrable T-periodic function, let 
( u *'')/eN be a sequence of non-vanishing piecewise 
constant T-periodic functions converging to u* in the 
distribution sense with ||u*''||l 1 ([o,t]) < |u|([0,T]) 

for every I. We define, for every n, u n j — , 

v n ,i ■ t h-> j' V\u n ,i\{r)dT and M n> i{t) = sg(u 

n.l ° 

v n ,i(t))e- Vn ^ A BWe v "'>W A . For every t, the matrix 
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Jq M/ i ;(r)dT tends to tM\ uniformly with respect to t 
in a compact set, as I tends to infinity. Hence the solu- 
tions of x = M l} i(t)x tend to the solutions of x = M^x, 
uniformly with respect to the time in a compact inter- 
val, as I tends to infinity. That concludes the proof of 
Proposition 2. 

2.4 Efficiency of the transfer 

Continuing with the notation of the last paragraph, for 
every non identically zero, A 2 ! r Afc | -periodic function u* , 
we define the efficiency of u* with respect to the transi- 
tion (j, k) as the real quantity: 

E {l - k) (u*) = 



j o ^-A fc i w * (r y (A ._ Afe)rdr 



Jo' 



\\j-X k \ 



u*(t)\At 



For every u, < E^' k \u) < 1. For every {j, k}, 
sup u i?^ ,fe '(u) = 1. The supremum is reached with a 
periodic sum of Dirac masses. 

An intuitive explanation of the efficiency could be the 
following: asymptotically, the L 1 norm of the control 
needed to induce the transition between levels j and 
k using periodic controls of the form u n is equal to 

n/(2\bj, k \E^\u*)). 

The system (2) being given, the design of an effective 
control law fulfilling the hypotheses of Theorem 1 is an 
important practical issue. To generate a transfer from 
level j to level k, one should choose a control u such that 
E^' k '(u) be as large as possible and E^ ll ' l2 \u) be zero 
(or arbitrarily close to zero) for every Zi , Z2 such that 
A; 1 — A; 2 G (Xj — Afc)Z. The algorithm we have described 
in [BCCS12] allows us to build u such that E^' k \u) > 
0.43, with E^ ll,l2 \u) arbitrarily small for every finite 
number of pairs {h,h} satisfying {h^h} ^ {j,k} and 
l-^i — 1 IAj — •M- Some other examples, including 
also examples of ineffective controls with zero efficiency, 
are studied in [BCCllb]. 

3 Infinite dimensional estimates 

We come back to the general case of Assumption 1. To 
ensure that the system (2) is well-posed, we consider 
only piecewise constant functions u* . The method of the 
proof is directly inspired by [BCCS12]: the original infi- 
nite dimensional system (2) is approached by a suitable 
Galerkin approximation, which allows us to apply the 
finite dimensional results of Section 2. 



ip E H '}Zj <N {4'j,'^)4>j € H and the compressions 
A( N ~) — ttn o A o tt n and B^ — ttn ° B o tt n . Note 
that A( N ^ and B^ are finite rank operators defined in 
an infinite dimensional space. With an obvious abuse of 
notation, we extend the finite dimensional propagator 
X^-j defined in Section 2.1 to the infinite dimensional 
space H. 

Let u* : R — > U, j, k be given as in the hypotheses of 
Theorem 1. Up to a reordering, we may assume j = 1 
and k = 2 without loss of generality. 



Define 



K 



2\b 12 \Eff(^)(u*)- 
Fix e > 0. Since and §1 belong to the domain of B, 
the sequences (6i,/);gn and (&2,/)zeN are in I 2 . Hence, 
there exists N in N such that ||7T2-B(1 — ttjv) || = ||(1 — 
ir N )Bir 2 \\ < e/K. 

Consider system (2) with control u n := u* /n in projec- 
tion on span(0i, . . . , 4>n). 



TTN 



dt 



T?»0i = (A (N ) + Un (i)B (Ar) )T^! 



+u n (t)7r JV B(l-7r JV )TJ , »^i. (15) 
From the variation of the constant, we get 

TTjvT?"^! 0)01 

+ Ju n (a)X^ ) (t, s)ir N B(l - irrfTpfads. (16) 
Project (16) on span(0i, (fo), and recall that 7Tjv7T2 = 

7T27TJV = 7T2 f° r N > 2: 

TraT?-^ =^^(4, 0)01 

+ Ju n (8)n i X$ ) (t,8)w N B(l-* N )??»<t> 1 d8. (17) 

Define, for every t, s in R, the bounded linear mapping 
fa, XM(t, a)] ':= n 2 oX^ } (t, s)-X$ ) (t, s)oir 2 . Equa- 
tion (17) reads 

- Ju n (s)X^ ) (t, s)n 2 B(l - TrjvJT^^ds 
+ Ju n {8)[n 2 ,X^ ) {t,a)]ir N B(l-w N )r^<f, 1 d8. (18) 



3.1 Galerkin approximation 

Let (A, B, U) satisfy Assumption 1. We define bjk = 
(<f>j , B<pk) for j, k in N. For every N in N, we define ttn 



3.2 Estimates of commutators 

Extend the definition of M f of Section 2.3 by M^fa 
for j > N and define the linear operator E N (t) : 
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X"£ (t,0)-e vl 11 (*) M 1 Since the commutator [tt 2 , Aft] = with T* = n/2, I = 4/|A 2 - Ai|, K = 2/\b 12 \ 



n 2 M ji — M^ir 2 vanishes, we have, for every t in R, 

|| faXfc (t, 0)] || = || [Tra, e^ 11 + iSfo (t)] || 
= \\[*2,E? N) (t)]\\<2\\E? N) {t)\\. 

Note also that, for every t in R, 

Wi^X^OMl = ||^ ) (0,t)[^ ) (t,0),7r 2 ]X ( ^ ) (0,t)|| 
<2||25fo(t)||. 

For every s, t in R, 

fa, *)] 

= ^ ^""j (i, 0)*^ (0, a) - X"^ (i, 0) (0, S )^ 2 

= X^ ) M)[tt 2 ,A^ ) (0, S )] + [tt 2 , OflX^O, a). 

Finally, we get, for every (s, i) in R, for every n, AT in N. 

^Xfcfaa)] <4\\E? N) (t)\\. (19) 



3.3 Proof of Theorem 1 

From (18) and (19), since ||tt 2 B(1 - w N )\\ < e/K, 

<e + 4\\E? N) (t)\\K\\n N B(l-ir N )\\. (20) 

From (13), sup t<t , n ( K ) l-E-^WII tends to zero as n tends 
to infinity. For n large enough, 

\\E? N) (nT*)\\< 



4K\\n N B(l -ttjv) 



and 1 - \(cj)2,T^(nT*)(t>i)\ < 2e. This completes the 
proof of Theorem 1. 

3.4 Bounded coupling 

In this Section, we consider the special case where B is 
bounded. The propagator u H> T u is defined on the set 
of locally integrable functions. 

Proposition 6 Let (A, B, R) satisfy Assumption 1. As- 
sume that B is bounded and that (1,2) is a non de- 



2tt 



and 



generate transition of(A,B). Define T — , x 

u* : t M> cos((A 2 — Xi)t). Then, for every n in N, there 
exists T* in (nT* - T, nT* + T) such that 



l-|(02,T^(r n *,O)0O| (l + C')\\n 2 B\\I 

1 + 2A-IIBII S n ' { ] 



a 



sup 

C?,fe)eA' 



sm 7T 



|A 2 -Ai| 



where A' = {(j,k) G {1, . . . , N} 2 such that {j, k} n 
{1,2} ^ and \Xj - X k \ < 3/2|A 2 - Ai| andb jk + 0}. 

Proof: Fix e > 0, and define = e/K. We apply the 
procedure of Section 3.1 and find N in N such that || (1 — 
ttn)Btt2\\ < n. By definition of N, for every n in N, for 

every * < T*, ||7r 2 T^i - ^Xfift^faW < e. 

A direct computation shows that, for every co > 3tt/T, 



£e^u*(t)dt 



sin(wT/2) 



uT 2 



oj 2 T 2 - 4tt 2 



<^<j. 

5 7T 



Hence, for every j, fc such that Xj — X k > 3/2 |A 2 — Ai|, 
(9) reads 

2\b jk \I 



K k (r)dr 



< 



n 



Estimate (4) becomes 



1 - |(0 2 ,X ( ^(r K *, Q)0Q| (l + C')\\ir 2 B\\I 
l + 2Jf||B|| - n 

or, by definition of N, 

i-|(0 2 ,r^(r*,o)0Q| <2p 1 (i + CQIk 2 g||J 

l + 2Jf||B|| " n 

Conclusion follows when e tends to zero. □ 

Proposition 6 is the translation in mathematical terms 
of the well-known fact that the difficulty of inducing 
a given transition between two eigenstates j and k of 
the free Hamiltonian is mainly due to eigenstates with 
energy close to Xj or Afe. In other words, in the case of 
bounded coupling, the convergence in Theorem 1 occurs 
independently of the high energy levels of A. These levels 
are rarely known precisely in practice. 

4 Example: Morse potential with truncated 
dipolar interaction 

The approximation procedure of Section 3 has been 
successfully applied to some classical models, see for 
instance Boussai'd, Caponigro & Chambrion [BCCllb] 
and [BCClla] for the rotation of a 2D molecule, the in- 
finite square potential well (with bounded coupling op- 
erator B) and a perturbation of the quantum harmonic 
oscillator (with unbounded potential B) . 
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We present below another type of quantum oscillator. Its 
main feature is that the spectrum of the free Hamilto- 
nian presents a continuous part and a discrete part. Very 
few controllability results are known for such systems 
with mixed spectrum (sec Mirrahimi [Mir09] for control- 
lability results based on a Lyapunov approach). 



K = 2/1(^0,^)1, A' - {(j,k) e {0,l,2,3} 2 |{j,fc} ± 
{0, 1}} and 



C" = sup 

U,k)eA' 



sin 7T 



Xj — Afe 
Ai — A 



4-1 Modeling 

We consider a diatomic molecule submitted to a time 
variable electric field with support in some bounded do- 
main. The potential energy of the molecule is modeled 
with the Morse potential (see Morse [Mor29]). With a 
suitable choice of units, the Schrodinger equation reads 



. dtp d 2 tp 



dt 



dx 2 



+ V(x)ip +u(t)W M (x)ip, x>0, 



where V : x n> a 2 [e~ 2 ^ x ~ x ^ - 2e~( x ~ Xe ^), and W M ■ 
x n> x if x < M and W M {x) = if x > M. The positive 
constants a and x e describe the physical properties of 
the molecule: the potential V reaches its minimum at 
x e , which can be considered as the equilibrium length 
of the molecule and a is related to the depth of the 
potential well. The constant M (large with respect to 
x e ) is the size of the region where the electric field is 
assumed to be active. With the notations of Section 1.3, 
H = L 2 ((0, oo), C), A = i(A - V) and B = -iW M - 

Define the integer N = [a — 1/2J. In the following, we 
assume that N > 4. The skew-adjoint operator A ad- 
mits a finite family simple eigenvalues (— i\ n )o<n<N , as- 
sociated with the eigenfuctions <ft n : x i-> e~ x x a P n (x) 
where P n is a polynomial function. For every n < N, 
A„ = — (a — n — 1/2) 2 . The spectrum of A contains also 
a continuous part — i[0, +oo). 



From Theorem 2.1, page 525, of [Kat95], for every 77 > 0, 
there exists a skew-adjoint operator A v such that A v 
admits a complete family of eigenvectors, A<f>i = A„ 



for every I < N and \\A — A tl \\ < rj. The pure point 
spectrum (— iA^) n€ N is everywhere dense in — i[0, +00). 

The system (A v , B, R) satisfies Assumption 1. For every 
locally integrable u, we denote with the propagator 

of = (A, + uB M )^. 

We choose an integer n > (1 + C")MI(1 + 2MK)/e 
and we define rj = e/ (nT* + T) . The transition (0, 1) of 
(A v , B) is non degenerate, Ai — A = — 2 (a + 1) and, for 
every I > 4, 



A + |(Ai - A )^ > A 4 - (\ + |(Ai - A ) 
>5a- 17 > 0. 



Proposition 6 applied to system ( A Vl B, R) gives the ex- 
istence of T* < nT*+T such that 1 - 1 (</>i , (T* ) </> ) | < 
e. Since \\A - A v \\ < 77, we have, for every t < nT* + T, 



r^(t,o)-T r p (t,o) 



< (nT* + T)n < e, 



and finally 1 - | (T*, O))0 O } I < 2e. 



4-2 Galerkin approximation 



5 Conclusion 



We aim to transfer the system from the first energy level 
to the second one. Since the family (4> n ) n <N is not an 
Hilbert basis of H , (A, B, R) does not satisfy Assump- 
tion 1 and we cannot apply directly Theorem 1. 

Note that (<po,B(f>i) tends to J R+ x<po(x)<j)i (x)dx / 
as M tends to infinity. Hence (4>o,B(f>i) ^ for M 
large enough. Moreover, for every i, to, A to — A; = (n — 
m)(-2a + 1 + (n + m)). 

From now on, we do the generic hypothesis that a ^ Q. 
In this case, A m — A; = A;/ — X m r implies n — m = 
n' — m' and n + m = n' + to', that is {n, to} = {n' , to'}. 
Moreover, for every {j,k}, | — A fc | £ (N\{0})(Ai - A ) 
implies {j, k} = {0, 1}. 

Fix e > 0. Inspired by Proposition 6, we define 
u* t 1 y cos((Ai - Xo)t),T* = tt/2, / = 4/|Ai - A |, 



The contribution of this paper has two elements. First, 
it proves the validity of the Rotational Wave Approxi- 
mation for infinite dimensional quantum systems with 
a pure point spectrum. The result can be partially ex- 
tended to systems with a mixed spectrum (including a 
discrete and a continuous part). The convergence results 
are accompanied by explicit error estimates which pro- 
vide explicit time estimates for the approximate con- 
trollability of bilinear infinite dimensional quantum sys- 
tems. Second, a notion of efficiency has been introduced. 
The efficiency is a measure of the i^-norm of the control 
needed to achieve a given transition between two eigen- 
states of the free Hamiltonian. 

The L^-norm of the control appears to play a central 
role in the analysis of bilinear quantum systems. This is 
slighty surprising, since one would rather expect that the 
L 2 -norm of the control field represents the energy given 
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to the system. It may be an artifact of the semi-classical 
model that disappears with a more realistic model that 
takes into account a quantized field. 

Among other topics, future analysis may concentrate on 
the design of time-efficient controls or on a systematic 
treatment of quantum systems when the free Hamilto- 
nian has a continuous spectrum. 
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